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Topological superconductors differ from topologically trivial ones for the presence of topologically 
protected zero-energy modes. To date, experimental evidence of topological superconductivity in 
nanostructures has been mainly obtained by measuring the zero-bias conductance peak via tunneling 
spectroscopy. Here, we propose an alternative and complementary experimental recipe to detect 
topological phase transitions in these systems. We show in fact that, for a finite-sized system with 
broken time-reversal symmetry, discontinuities in the Josephson current-phase relation correspond 
to the presence of zero-energy modes and to a change in the fermion parity of the groundstate. 

Such discontinuities can be experimentally revealed by a characteristic temperature dependence of 
the current, and can be related to a finite anomalous current at zero phase in systems with broken 
phase-inversion symmetry. 


Introduction The recent discovery of topological ma¬ 
terials deeply impacted condensed matter research[T]. 
These materials exhibit a number of exceptional prop¬ 
erties which are related to the presence of topologi¬ 
cally protected states localized at their edges. In topo¬ 
logical superconductors (TS), for instance, Majorana 
edge states[2] are characterized by a distinctive non- 
abelian statistics, which makes them ideal candidates for 
fault-tolerant quantum computation [3]. Theoretically, a 
TS can be realized at the interface between a conven¬ 
tional superconductor and a topological insulator [1], or 
in semiconductor-superconductor heterostructures with 
spin-orbit coupling (SOC) in magnetic field [5]. 

However, revealing signatures of topological non-trivial 
phases in TS’s is not straightforward [6]. This is mainly 
because, unlike more conventional continuous phase tran¬ 
sitions, topological phase transitions [7] do not break any 
symmetry nor exhibit any critical behavior, but are in¬ 
stead identified by a change of the corresponding topolog¬ 
ical invariant [8] and of the edge properties in continuous 
systemsj^. To date, experimental evidence of non-trivial 
superconductivity comprises the measure of the 47r- 
periodic Josephson current in TS rings [TU|. or the zero- 
bias conductance peak via tunneling SDectroscoDv|llj or 
spatially-resolved spectroscopic imaging [I^. 

In this work we propose a direct method to probe vari¬ 
ations of the topological invariant, i.e., the fermion par¬ 
ity of the groundstateflB]. in Josephson junctions [HI ITS] . 
It is known that the Josephson current-phase relation 
(CPR) may exhibit discontinuities at low temperatures 
corresponding to Andreev level crossings [H], in systems 
as different as quantum dots|T7], nanowires [IB]. Joseph¬ 
son junction arrays [TB], and Weyl semimetals [BO]. Here, 
we show that in hnite-sized non-interacting s-wave TS’s 
with broken time-reversal symmetry, such discontinuities 
are related to the presence of zero-energy modes and to 
variations of the fermion parity, which define the topo¬ 
logical phase transitions of this system. These discon¬ 
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FIG. 1. A TS (a) realized by a quantum nanostructure 
embedded into a superconducting ring with Rashba SOC 
oc (cr X p)^ along the « axis. The nanostructure (b) of length 
L, width IT, and proximity length Lp {Ln = L — 2Lp), con¬ 
nected to two superconducting leads, can be either a quantum 
dot (L = W = 1 lattice site), wire (L > W = 1), or planar 
well {L,W> 1). 


tinuities can be experimentally revealed by a character¬ 
istic temperature dependence and are moreover related, 
in systems with broken phase-inversion symmetry, to an 
anomalous current at zero phase. 

The model A TS can be realized by a nanostructure 
embedded into a superconducting ring[2T] with SOC and 
magnetic field [5] (see Fig. [^. Here, we consider the 
case of hnite-sized systems, e.g., a quantum dot, wire, 
or planar well (see Appendix 0, described by a tight- 
binding Bogoliubov-de Gennes (BdG) Hamiltonian [B5] 
with linear dimensions smaller than the coherence length 
^ = hvp/A, where vp is the Fermi velocity and A the su¬ 
perconducting gap. This system is the zero-dimensional 
(OD) limit of a one-dimensional (ID) Majorana chain[2], 
in the sense that its Hamiltonian does not depend on 
any momentum-like continuous parameter. The Andreev 
spectrum of this system is a discrete set of particle-hole 
energy levels Ei (tp) which depend on the gauge-invariant 
phase difference ip of the superconducting order param¬ 
eter between the two leads, induced by a magnetic flux 
$. If L <C where L is the distance between the leads, 
the CPR is given bv|24j 
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FIG. 2. Topological phase space (a) of a finite-sized TS 
as a function of </?, 6 and A. Topological transitions between 
states with different fermion parity are possible for |A| < 1, 
where Vtp = ±1 respectively for \ip — 6\^ arccos (—A). Zero- 
energy modes occur at cos(i/p — 9) — —A (solid line). CPR in 
units of /o = {e/h)A, LE Andreev level, and Pfaffian T^jC 
for a quantum wire (b) and a planar well (c) with magnetic 
field b = by perpendicular to the SOC and current directions. 
Different magnetic field directions give qualitatively similar 
results. For |A| < 1 (continuous lines), CPR discontinuities 
correspond to zero-energy modes (J>* = 0 ) to variations 
of the fermion parity Vip = sgn at = 0 ± arccos(—A). 
The trivial Tp > 0 and non-trivial Ty, < 0 branches of the 
CPR correspond respectively to A > 1 (dashed lines) and 
A < — 1 (dotted lines). In a planar well with by ^ 0 and 
a 7 ^ 0 (c), the CPR and Pfaffian are no longer symmetric 
under phase-inversion —>■ —ip (9 7 ^ nn in this case), and an 
anomalous current /(O) 7 ^ 0 is present for A = —0.9 and —1.5. 


where f{x) = l/(e^ -I- 1) is the Fermi-Dirac distri¬ 
bution. Note that the low-energy Andreev spectrum 
does not depend on the superconducting ring length, 
in the limit of short iunctions|25j (see Appendix g. 
If the lowest energy (LE) levels close the gap with lin¬ 
ear phase-dispersion, the CPR exhibits a discontinu¬ 
ity at zero temperature T = 0. In fact, in this case 
the Fermi-Dirac distribution in Eq. Q converges to a 
step function, and thus the only contribution to the 
current is given by energy levels E < 0, i.e., I{f) = 
This mandates a discontinuous 
drop = -{e/h)Y,E,(v-)=o \dvE 3 {.‘P *)\at any gap¬ 

less point <p* where zero-energy levels Ej{if*) — 0 have 
linear phase-dispersion d^Ej{(p*) 7^ 0. Since Andreev 
levels are continuously differentiable in finite-sized sys¬ 
tems (see Appendix g and Refs. |26]) these are the only 
points where the CPR can be discontinuous. Hence, 
discontinuities at zero temperature correspond to zero- 
energy modes closing the particle-hole gap. In general, 
the converse is not true, i.e., the gap may close with¬ 
out any discontinuity if dipEj{ip*) = 0. Hereafter we 
will show that such discontinuities correspond, if time- 
reversal symmetry is broken, to fermion parity transitions 
and thus can be related to a topological phase transition 
of the TS. 

Superconductors exhibit particle-hole symmetry, i.e., 
their BdG Hamiltonian ’H(</3) is invariant under the an¬ 
tiunitary transformation S = TxK with = 1, where 
K is the complex conjugate operator and Tx the Pauli 


matrix in particle-hole space. For a finite magnetic field 
b 7^ 0, time-reversal and chiral symmetries are broken, 
and hence the system is in the Altland-Zirnbauer [27] 
symmetry class D. This class is characterized both in 
ID (Majorana chain, continuous spectrum) and in OD 
(finite-sized system with discrete energy spectrum) by 
the Z2 topological invariant accordingly to the peri¬ 
odic table of topological phases [27]. Analogously to 
the ID case (Majorana chain), the topological invari¬ 
ant in a OD system is defined following Ref. [25| as the 
fermion parity of the grounds!at ejH |T3] Vy, = sgn 
where Ey, = ^i\H{p>)iTx\ is the Pfaffian [2^ of the ma¬ 
trix 'H{ip)iTx- The fermion parity labels the topologi¬ 
cal inequivalent phases, i.e., trivial V^p = 1 > 0) 

and non-trivial = — 1 {E^p < 0). Moreover, since 
E'^ = Aet[T-L{ip)iTx] = det[H(v5)] = ]\iEi{ip), the con¬ 
dition E^* = 0 corresponds to gapless points ip* where 
zero-energy modes occur. The topological invariant in 
OD (finite-sized system) and in ID (continuous limit) are 
closely related: in the limit L —>■ 00 in fact, the Majo¬ 
rana number |2] coincides with the fermion parity of the 
groundstate, i.e., M. = senipf ['Hzr.,.]l [OlE5] . 

Fermion parity phase-dependence The non-trivial 
phase of a Majorana chain (ID, continuous spectrum) 
requires an open particle-hole gap, which can be real¬ 
ized only with SOCjSj and for specific magnetic field 
directions|30j. Contrarily, in OD finite-sized systems, the 
topological invariant Vip is well defined even in absence 
of SOC and for any magnetic field direction, as long as 
E^p 7^ 0, and depends explicitly on the phase ip. Changes 
in the fermion parity Vp, define the topological phase 
transitions in this system. Independently from the details 
of the tight-binding Hamiltonian, the Pfaffian can be ex¬ 
panded as a Fourier series in the phase p with coefficients 
On oc A^". If the superconducting gap is smaller than the 
bandwidth of the nanostructure (e.g., in conventional su¬ 
perconductors), higher harmonics become negligible and 
one obtains at the first order (see Appendix [b|) 

E^^C [cos((^ - 0) + A], (2) 


where A = {Ep Eji^C) (the sum does 
not depend on the choice of the angle y>'), 2(7 = 
i/(J'o - -^77)^ + (-^V2 “ -^-72)^, and tan0 = (E^j^ - 
E_„/ 2 ) I {Eq — E-rr)- These parameters depend on the 
Hamiltonian, e.g., on the hopping parameter t, magnetic 
field b, chemical potential /r, SOC a, and superconduct¬ 
ing gap A, but not on the gauge-invariant phase p. The 
topological phase space is thus completely characterized 
by the phase p and the parameter A = A(t, a, p, b. A), as 
shown in Fig. 2(a) If |A| > 1 the system is either in the 
trivial Ep > 0 for A > 1 or non-trivial phase Ep < 0 for 
A < — 1, with the exception of a single gapless point at 
p* — 9 TT or 0 respectively for A = ±1. If |A| < I in¬ 
stead, topological transitions occur at the gapless points 
p’^ ^ 9 ± arccos(—A). Hence, the closing of the particle- 
hole gap Ep = 0 defines the boundaries between trivial 

























3 


and non-trivial phases. 

Furthermore, if SOC or magnetic field vanish, the An¬ 
dreev spectrum, CPR, and fermion parity are invariant 
under phase-inversion, which mandates B = titt with n 
integer in Eq. (§■ This invariance is related to the mag¬ 
netic mirror symmetry [3T] = K, where M^z = ^o^y 

is the spin mirror reflection across the xz plane and 
0 = —layK the time-reversal operator. If a = 0 or b = 0 
in fact, the quantization axis can be arbitrarily chosen 
such that the only complex terms in the BdG Hamilto¬ 
nian are those in the phase ip (see Appendix . There¬ 
fore 'H{ip)* = T-Li—ip), and consequently I{ip) = —I{—(p), 
Vip = V-ip, and 9 = mr. However, if the magnetic mirror 
symmetry is broken, i.e., if 'H{(p)* ^ T-Li—ip), the CPR 
and the fermion parity may be no longer symmetric un¬ 
der phase-inversion, i.e., I{(p) ^ —I{—ip), Vp ^ 'P-p, 
and 9 ^ mr. In this case, the magnetic mirror symme¬ 
try corresponds to the inversion of the Pfafhan phase- 
shift 9 —)■ —9 (see Appendix and the fermion par¬ 
ity is still invariant under the more general transforma¬ 
tion ip ^ 29 - ip, i.e., Vp = V 2 e-p [cf. Eq. 0]. A 
Pfafhan phase-shift 9 ^ mr \s thus a signature of the 
broken phase-inversion and magnetic mirror symmetries, 
and can result in an anomalous current p2l434) at zero 
phase (see below). 

CPR discontinuities In a superconductor with bro¬ 
ken time-reversal symmetry, any zero-energy mode is at 
least doubly-degenerate due to particle-hole symmetry. 
Hence the Pfafhan can be expanded near any gapless 
point as iFp oc {(p — ip*)‘^ where the order d is half the 
total multiplicity 2d — ruj , with ruj the multiplici¬ 
ties of zero-energy modes Ej{ip) (x {ip — Thus, 

if the Pfafhan hrst derivative = dpEp is non-zero 
at the gapless point ip*, there exist only two doubly- 
degenerate LE levels E±{ip*) = 0 with linear phase- 
dispersion E±{ip) oc {ip — ip*). Therefore one obtains 
dpE±{ip*) = ±E{p./x^-, where = UE,>oEiiP*) > 
0 (see Appendix [ 0 . Hence, the LE level contribution 
changes its sign passing through the gapless point, and 
the total current exhibits a discontinuous drop given by 

= (3) 

n xv* 

where \ = C'|sin((/?* — 0)| = Cy/l — [cf. Eq. 0]. 
Equation 0 relates CPR discontinuities at zero tem¬ 
perature with the variations of the fermion parity Vp = 
sgaEp in superconductors with broken time-reversal 
symmetry, and is valid for F'p, ^ 0. If |A| < 1 in fact, 
the Pfafhan changes its sign at ip"^ = 9 ± arccos(—A), 
where ^ 0 according to Eq. Q. Here, the CPR 
has two discontinuities A/((^ 5 i) 0 which correspond 

to transitions between topological phases with even and 
odd fermion parity, where zero-energy modes appear. If 
|A| > 1 instead, no zero-energy mode nor CPR discon¬ 
tinuity occur. In the limit cases A = ±1, the Pfafhan 


vanishes at the gapless point without any sign-change 
{Vp, = 0), while the CPR may exhibit at most one 
discontinuity, corresponding to a zero-energy mode with 
m± = 1. Hence, the presence of two distinct discon¬ 
tinuities in the CPR of a TS with broken time-reversal 
symmetry dehnes the boundaries between inequivalent 
topological phases. Moreover, these discontinuities are a 
direct signature of zero-energy modes, and coincide with 
a sign-change of the LE level contribution to the current. 
These zero-energy modes signal the topological transi¬ 
tion at ip* between phases with different fermion par¬ 
ity, and are described locally as a linear superposition 
of particle and hole states or, equivalently, of two or¬ 
thogonal Majorana states (see Appendix |0. The results 
presented here hold for any discrete OD Hamiltonian in 
the Altland-Zirnbauer class D (particle-hole symmetry 
and broken time-reversal symmetry). Note that, if time- 
reversal symmetry is unbroken (b = 0 and 9 = mr), the 
total multiplicity is 2d > 4 due to spin and particle-hole 
degeneracy, and therefore = 0 at any gapless point. 
Hence, no fermion parity transition occurs (|A| > 1) 
and the CPR may exhibit at most one discontinuity at 
ip* = mr if |A| = 1, e.g., in point contact iunctions |24). 
Indeed, a time-reversal invariant s-wave superconductor 
is topologically trivial. Since CPR discontinuities corre¬ 
spond to an abrupt jump to the LE state, they can be 
measured at the equilibrium (DC Josephson current) and 
are not affected by quasiparticle poisoning|35). 

Numerical results Figures |2(b)| and |2(c)| show the 
CPR at zero temperature, the LE Andreev level, and the 
Pfafhan of a quantum wire {W = 1) and a planar well 
{W > 1) as a function of the phase and magnetic held 
b = by perpendicular to the SOC [oc (cr x p)^] and cur¬ 
rent directions, calculated directly from the BdC Hamil¬ 
tonian ^{{ip) (see Appendix [A| for details). Eor |A| < 1, 
the gap closes with linear phase-dispersion and thus the 
CPR has two discontinuities at ip'^. In a planar well 
with by ^ 0 and a ^ 0 [Pig. |2(c) , the CPR and the 
Pfafhan are no longer symmetric under phase-inversion 
ip —>■ —ip {9 ^ mr in this case). Note that fermion par¬ 
ity transitions and CPR discontinuities are present for 
any magnetic held direction. In a quantum wire with 
L —> 00 (ID, continuous spectrum), we verihed numer¬ 
ically that Eqs. 0 and 0 reproduce the well-known 
results of Refs. El In particular, the trivial {Ni = 1) and 
non-trivial {A4 = — 1 ) phases correspond respectively to 
A > 1 and A = — 1 in Eq. 0 . We have also verihed 
that Majorana bound states localize at the wire edges, 
corresponding to a discontinuity in the CPR at ip* = ir 
it h ± y. 

Anomalous current When the magnetic mirror sym¬ 
metry is broken {9 ^ mr), the current may be no longer 
symmetric under phase-inversion, which can result in a 
hnite anomalous current [5^34j at ip = 0, as shown in 


Eig. 2(c) This can be realized, e.g., in planar quantum 
wells with SOC and magnetic held, where a hnite anoma- 
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lous current has been related to the presence of chiral 
edge states[33] or to a SOC-induced Lorentz force [34], 
In these systems, a topological phase transition can be 
revealed by a discontinuity of the current at zero phase 
AI(ip = 0, v) with respect to the parameter v (magnetic 
field, chemical potential, or SOC) which drives the sys¬ 
tem through the topological transition. In this case one 
can always find a value v = v* such that the gap closes 
at if* = 0, i.e., \(v*) = — cosO [cf. Eq. ([^ and Fig. 2(a)| . 
At this gapless point the fermion parity changes and the 
current at (/? = 0 has a discontinuous drop with respect 
to the parameter v given by 

AI{ip* = Q.v = V*) =\sin0\. (4) 

« Xo 

The value v* corresponds to a crossover between the triv¬ 
ial > 0 (A > 1) and non-trivial JA< 0 (A < —1) 


branches of the CPR, as shown in Fig. 2(c) This discon¬ 


tinuity mandates a finite anomalous current 7(0,;/) 7 ^ 0 
near v = v* either in the trivial or non-trivial phases, or 
in both. Numerical calculations indicate that it is non¬ 
zero only in the non-trivial phase [cf. Fig. 2(c)| . Hence, a 
discontinuity of the anomalous current at zero phase with 
respect to any system parameter (e.g., magnetic field) is 
also a signature of a topological phase transition. 

Experimental proposal At finite temperatures, CPR 
discontinuities are smoothed out by the thermal spread¬ 
ing of the Fermi-Dirac distribution. For T = 

Sd/ks, where 6d is the gap between the first and sec¬ 
ond Andreev levels at (p*, the current can be expanded 
as a sum of two contributions Ihe{^) and hei'p) coming 
respectively from higher energy levels {E > 8d) and from 
the LE level. The latter contribution depends strongly 
on the temperature and can be obtained from Eqs. Q 
and ([^ (see Appendix [Pj) , which yield 


I{ip)^Ihe{p)+ ^^^!f ^ tanh 


h AI{(p*){ip - ip*) 
e 4:kBT 


( 5 ) 


The current phase-derivative diverges as « 

—H/{8e)AI{ip*)^/{ksT) for T —0. Hence, the scal¬ 
ing factor Sip* = —dipI{(p*)kBT is a direct measure of 
the discontinuous drop, since AI{(p*) « —*yj8ejK)s^. 
Since the separation between energy levels increases as 
the system size decreases, the temperature T^, below 
which CPR discontinuities are measurable, is maximized 
for small linear dimensions. For parameters considered 
in Figs. and we obtain Td « 0.2 Tc. 

In light of this, we propose to measure the low- 
temperature CPR[M1 13S] through a nanostructure with 
broken time-reversal symmetry. Figures |3(a )] and |3(b)| 
show the CPR of a quantum wire (L = 200 sites) at 
T > 0 and its phase-derivative as a function of the mag¬ 


netic field. Besides, Fig. 3(c) shows the effect of temper¬ 


ature on current discontinuities, which is described by 
Eq. (§. Spikes in the phase-derivative which exhibit the 



FIG. 3. Fermion parity Vip, CPR (a), and phase-derivative 
dipl{ip) of the current (b) at T = 0.02 Tc of a quantum wire 
as a function of the magnetic held along the y axis. Spikes in 
the phase-derivative correspond to topological phase transi¬ 
tions. Josephson current (c) near the gapless point and phase- 
derivative 7 ((^*) (inset). Dotted lines correspond to Eq. ([^, 
continuous lines to numerical calculations. The two minima 
(b) at i/p = 0 {by/ A ~ 0.5 and « 2 respectively) are constant 
in temperature, i.e., do not diverge for T —0, and indeed do 
not correspond to any topological phase transition. 


characteristic temperature scaling of Fig. 3(c) identify 
the boundaries between inequivalent topological phases. 
These effects should be measurable in clean Josephson 
weak-links of the order of 50-150 nm in InAs or InSb 
nanowires {L = 100-300 lattice sites) proximized by a 
conventional superconductor [3T] (e.g., Nb or Al). Note 
that the Josephson critical current has been recently mea¬ 
sured in an InAs nanowire in the relevant regime [22]. at 
very low temperatures ~ 10 mK comparable with the 
temperatures considered here. 


Conclusions We have shown that discontinuities in 
the Josephson CPR correspond, in TS’s with broken 
time-reversal symmetry, to topological phase transitions 
between states with different fermion parity, where zero- 
energy modes occur. These current discontinuities are 
not affected by quasiparticle poisoning and can be re¬ 
vealed by spikes in the current phase-derivative and by 
their characteristic temperature dependence. Moreover, 
in systems with broken phase-inversion symmetry, topo¬ 
logical phase transitions correspond to discontinuities of 
the anomalous current at zero phase. Such features in the 
CPR provide an experimental tool to probe the fermion 
parity and to resolve the topological phase space of TS’s. 
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Appendix A: Bogololiubov-de Gennes Hamiltonian 


In general, a ballistic Josephson iunctionfli] realized 
by a conventional superconductor in a magnetic field can 
be described by a tight-binding BdG Hamiltonian ^5] 
with broken time-reversal symmetry and with an s-wave 
spin-singlet superconducting pairing. In order to perform 
actual calculations, we consider specihcally the following 
Hamiltonian 




bi • cr - /rcTo 




-(bi • (T - ^CTo)’’ 






tcTQ (Ty 0 

0 -(tao+iaijay) ^ 


0 -{t(7Q+iaija^) 


-I- h.c., 
+h.c., (AI) 


where ’®'i = (c^^, c^, is the Nambu spinor, 

the Pauli matrices, fi the chemical potential, t the hop¬ 
ping parameter, b^ = the Zeeman magnetic 

field, aij the Rashba SOC, and Ai{ip) the s-wave super¬ 
conducting order parameter. We assume uniform mag¬ 
netic field bi = b and SOC = a > 0 inside the 
nanostructure (of length L and width W in units of lat¬ 
tice sites); b^ = aij = 0 in the superconducting leads; 
Ai(ip) = 0 in the normal region of the nanostructure (of 
length Lff ); Ai(vj) = A and Ae®'^ respectively in the 
two superconducting leads (of length Ls) and in the two 
contiguous proximized regions (of length Lp). The coor¬ 
dinate system is chosen such that the quantization axis 
is along the z direction, the supercurrent flows along the 
X direction, and the Rashba SOC acts along the z di¬ 
rection, being oc (cr x p)^. The BdC Hamiltonian AI 
can describe a quantum dot {L = W = 1 lattice site), a 
wire {L > W = 1), or a planar well {L,W > 1). More¬ 
over, we assume that, in general, the nanostructure is 
smaller than the superconducting leads, i.e., L Lg. 
Figure [AI| shows the energy spectrum of a quantum wire 
{W = I) with Ln = 100, Lp = 50, a = O.H, A = O.OOSt, 
^ = 2t, and for for b = 0 and b/A = (0.3,0.3,0.3) 
as a function of the phase (p and of the length of the 
superconducting leads Lg. The number of energy lev¬ 
els increases as the length of the superconducting leads 
Lg increases, but the low-energy Andreev spectrum (and 
in particular the LE level) remains unchanged (see the 


rightmost panels of Fig. All. In this case in fact, the 


degrees of freedom of the superconducting leads may be 
integrated out, and consequently the low-energy Andreev 
spectrum can be described by an effective tight-binding 
Hamiltonian in a finite-dimensional Hilbert space where 
the superconducting order parameter within the leads is 
replaced by an effective order parameter at the boundary 
of the nanostructure [5^. Hence, the low-energy Andreev 
spectrum does not depend on the ring length. There¬ 
fore in our numerical analysis we consider the limit of 



FIG. AI. Andreev spectrum of a quantum wire with 
Ln = 100, Lp — 50, a — O.lt, A = O.OOSt, p — 2t, for 
b = 0 and b/A = (0.3,0.3,0.3) respectively with |A| > 1 
(above) and |A| < 1 (below), as a function of the phase p 
with Lg = 50,100, 200, and as a function of the length Lg 
of the superconducting leads (rightmost panels) a.t p = n 
(|A| > 1, above) and at the gapless point p* (|A| < 1 below) 
respectively. The LE level (in red) and the gap 5d between 
the LE and the second lowest energy level do not depend on 
the length Lg. 


infinite ring length. The analytical results obtained in 
this work are largely independent on the details of the 
Hamiltonian |Al| and in particular do not depend on the 
boundary conditions (open or close) or on the presence 
of small perturbations (e.g., disorder). 


The Andreev spectrum of a confined electron system 
is a discrete set of energy levels Ei(p) which depends on 
the gauge-invariant phase difference p. Moreover, the 
Andreev levels of the Hamiltonian |A1| are continuously 
differentiable functions of the gauge-invariant phase dif¬ 
ference p. In fact, the energy levels are the roots of 
the characteristic polynomial of the effective low-energy 
Hamiltonian, whose coefficients depend smoothly on the 
gauge-invariant phase difference p. As shown in Refs. 1^51 
any polynomial whose roots are real and whose coeffi¬ 
cients are smooth functions of a parameter admits roots 
which are continuously differentiable in that parameter. 
Therefore the Andreev levels of the Hamiltonian ED are 
continuously differentiable in the phase p, for a suitable 
ordering which removes level-crossing discontinuities. 


Figures!^ and 1^ have been obtained calculating the An¬ 
dreev spectrum, the Josephson current, and the Pfaffian 
of the system described by the Hamiltonian m and as¬ 
suming open boundary conditions. In particular for the 
quantum wire [Figs. 2(b) ^and El we use Ln = 100, 
Lp = 50, W = 1, a = O.lt, A = 0.003t, p = 2t, and 
hy/lS. = 0, 0.2, 0.3, 0.4, 0.5, and 0.6, which correspond 
respectively to A = 1.1, 0.7, 0.2, -0.3, -0.7, and -1.1. For 
the planar quantum well [Fig. 2(c)| we use instead and 


Ln = 20, Lp = 10, IF = 4, a = 0.3t, A = 0.003t, p = 2t, 
and by/A = 14.0, 14.4, 14.7, 15.0, 15.2, and 15.6, which 
correspond respectively to A = 1.2, 0.4, -0.1, -0.6, -0.9, 
and -1.5. Note that the lattice parameters of InAs and 
InSb nanowires are « 4.5 A, which corresponds to a total 
length of « 100 nm for L = 2Lp + Ln = 200 lattice sites. 
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Appendix B: Phase dependence of the Pfafflan 

The invariant is defined in terms of the Pfaffian[29] 
of the BdG Hamiltonian in Majorana representation, 
which is the real number defined as J-^ = pf \}i{{p)iTx\- 
Since the Hamiltonian is periodic in the order parameter 
phase (/?, the Pfaffian can be expanded as a Fourier series 
as 


= Ao-^'^AnCos{nip-Bn) = ^ 




(HI) 


n— — oo 


where and Cn are respectively real and complex fac¬ 
tors, with Co = ^0 and C„ = respectively for 

n ^ 0. The tight-binding BdG Hamiltonian of an s-wave 
TS can be written in Majorana representation as 


A? 




, (A^)* of the Hamiltonian ap¬ 
pears at most once in Eq. ( |B4[ ) according to the proper¬ 
ties of the Pfaffian. Since the Hamiltonian |B2] is invariant 
up to a global change of the order parameter phase (i.e., 
Ai —>■ A^e®^) and we assumed A^ = A and A^ = Ae®*^ in 
the two superconducting leads respectively, the Pfaffian 
can depend only on the gauge-invariant phase difference 
Lp between the leads (or on its multiples mp). Hence, each 
term of the Pfaffian expansion Eq. (B4) has to be pro¬ 


portional to with n having its maximum value 

n = N in the case that Pk = Qk = ^ for any k. Therefore 


the Fourier expansion Eq. (Bl| has a finite number of 
terms, and is given by 


N 




= Ao + '^An cos {np - 9n), 


(B5) 


n—1 


n{p)iT^= - 


layAi 


ICJyAM 

—iayA\ 


Ha 


-m 




, (B2) 


where Hq does not depend on the order parameter, and 
N is the total number of lattice sites. Note that 'H{ip)iTx 
is indeed an antisymmetric matrix, since the terms layAi 
are antisymmetric. The Pfaffian of an antisymmetric ma¬ 
trix A of order 2n is defined as the polynomial of matrix 
entries given by 


Pf(^) = (B3) 

TT 

where the sum is over all permutations in the form 


1 2 3 
ii ji *2 


2n 

jn 


and such as ik < jk and ik < ik+i for any index k, and 
where sgn(7r) is the sign of the permutation (see, e.g.. 
Ref. [2^ . The Pfaffian is related to the determinant of 
antisymmetric matrices by the relation pf(H)^ = det(A). 
For the properties of permutations, each matrix element 
of 'H{p)iTx can appear only once in any of the product in 
Eq. (|B3[). Moreover, for the additional restriction ik < jk 


on the permutations in Eq. (B3), one has to consider only 
products of matrix elements above the diagonal. As a 
consequence each term of the Pfaffian can either depend 
linearly or not depend on the matrix elements A^ and 
A*. Therefore, one can expand the Pfaffian Ty, as 


_p — 


Pk,gk=o 


Af ... A^« (Af ... A?f )* , (B4) 


where Vql.-.qjj are complex prefactors which do not de¬ 
pend on the order parameter. Note that each term 


with An (X A^". If one defines 2C = max — minFq, 
and A = (maxJ^^ -|- min Jf^)/(2C'), one obtains 


= c 


N 


oo + a„ cos {np — On) + A 


n—1 


(B6) 


where ag = AgfC — A, a„ = An/C for n > 1, and 
/v = J-/C — A. It follows that \ ftp\ <l and thus Tq, = 0 


only if |A| < 1. Since the expansion in Eq. (B6) has a 
finite number or nonzero Fourier coefficients, the Pfaf¬ 
fian Tq, is an analytic and periodic function of the phase 
ip, and therefore it has a finite number of zeros in the 
interval [—7r,7r]. The same argument can be applied to 
the derivatives of any order, i.e., the Pfaffian has a finite 
number of local maxima, minima, and of saddle points 
(if any), which are indeed isolated points in the interval 

[-7r,7r]. 

If the superconducting gap A is small with respect to 
the relevant energy scales of the system, i.e., if A <C t, 
Eq. (B6| can be further simplified. In this case, since 
a„ oc A^", higher order Fourier coefficients can be safely 
neglected, and one obtains Tq, ~ Ag + Ai cos {p — Oi), 
which leads to Eq. ([^, i.e., TyjK, C [cos ipp — 0) + A]. In 
this approximation, it can be easily shown using trigono¬ 
metric identities that the coefficients A, 0, and C are 
given by 


.. ‘A p 4- J~ 

2 C ^ 


tan0 = 


-A ^/2 - 

-An — 


72 


2 C = J{Fg - + [T.U - 


(B7) 


where p' is any angle. The phase-shift of the Pfaffian 0 


in Eq. (B7) is related to the phase-inversion symmetry 
in the system. If the SOG or the magnetic field van¬ 
ish, the system is invariant under phase-inversion sym¬ 
metry, i.e., Ei{p) = Ei{—p) and = E-q,, which man¬ 


dates 6 = rni in Eq. (B7). This invariance is related to 
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the magnetic mirror symmetry ODerator|31j = K, 

where = lay is the spin mirror reflection across the 
xz plane and 0 = —layK the time-reversal operator. 
In general, the only terms in the BdG Hamiltonian |A1| 
which are affected by complex conjugation are the order 
parameter term, the Rashba SOC oc tax, and the term in 
the magnetic field component by. However, if a = 0, the 
only complex terms of the Hamiltonian |Al| are those in 
the phase ip and in the magnetic field component by. In 
this case, since the choice of the quantization axis is arbi¬ 
trary, one can always choose the z axis in the direction of 
the applied magnetic field. On the other hand, if b = 0, 
the only complex terms of the Hamiltonian are those in 
the phase and the Rashba SOC term oc tax, which trans¬ 
forms as (tax)* = —tax- In this case, a rotation of the 
quantization axis around the axis y is sufficient to re¬ 
store the sign of the SOC term. Therefore if a = 0 or 
b = 0 one has = 77(—i^)* and, consequently, the 

energy spectrum and the Pfaffian are symmetric under 
phase-inversion, i.e., Ei{ip) = Ei(—(p) and Ty, = T-y,. 
Therefore if a = 0 or b = 0 one has J-y, = J--y and 
Eq. (B7) further simplifies to 


ing the crossover between open and close boundaries via 
a weak link at the edge of the system. 


Appendix C: CPR at zero temperature 

In systems with particle-hole symmetry, the Andreev 
spectrum is particle-hole symmetric, and one has 

det [n{p>)\ = \{E,{yt) = n E,{v?. (Cl) 

Ei Ei>^ 

We note that the product on the right side of this 
equation is positive since the total number of nega¬ 
tive Andreev levels is even, due to the spin degree of 
freedom. Moreover, any zero-energy mode is at least 
doubly-degenerate, and can be in general expanded as 
Ej{ip) oc {ip — ip* near the gapless point ip*, where 
rrij > 1 is the multiplicity of the jth LE level at zero 
energy, which gives 

det [H(<^)] cx n(‘^ - = (T’ - , (C2) 

3 


A = 




■/2 




72 


C C 
2C=\Eq-E^\. 


9 = 


0 Eo> Et, 

TT E n ^ E rr 


(B8) 


On the other hand, if the magnetic mirror symmetry 
is broken, i.e., if Hiip) ^ Hi—ip)*, the energy spec¬ 
trum and the Pfaffian may be no longer symmetric un¬ 
der phase-inversion. In this case, the Fourier coefficients 


in Eq. (Bl) transforms under complex conjugation as 
Cn —>■ C* which mandates —>■ —Therefore, the 
magnetic mirror symmetry K corresponds to the inver¬ 
sion of the phase-shift of the Pfaffian 6 —)■ —9 in Eq. (§ 
and, as a consequence, the Pfaffian is no longer invariant 
under phase-inversion ip ^ —ip alone, but is still invari¬ 
ant under the more general transformation ip —3 29 — ip, 
I.e., E'_p E‘2.33 — y. 

There is a close relationship between the topological 
invariant in OD (finite-sized system), i.e., the fermion 
parity of the groundstate Vy = sgnJT^, and the topo¬ 
logical invariant in ID (continuous limit), i.e., the Ma- 
jorana number[2] AI. In a system with closed boundary 
conditions, any gapped phase is characterized by the ID 
topological invariant 


M. = sgnjpf [ 771 X 3 ,]} 

= sgnjpf [U{Q)iTx] pf [U{tt)iTx\}, 


(B9) 


where 77(fc) is the Fourier transform of the Hamiltonian 
77, as shown in Refs. |2] and [131 From Eq. (|B9|) follows im¬ 


mediately that the Majorana number M. coincides with 
the fermion parity in a system with closed boundary con¬ 
ditions. In a system with open boundary conditions in¬ 
stead, a similar relationship can be obtained by consider¬ 


where the index j spans only the modes which have 
zero-energy ai, ip = ip*, i.e., Ej{p*) = 0. Furthermore, 
for the relationship between the determinant and the 
Pfaffian one obtains that det[77((p)] = det [77 (<p) 1 X 3 ,] = 
pf[ 77 ((/ 3 )iX 3 ,]^, which yields 

Eyix{p-p*)-^^^^^, (C3) 

where Ey = \3i\H{p)iTx\ is the Pfaffian in Majorana rep¬ 
resentation. Therefore, the lowest order d of the non-zero 
terms in the Pfaffian expansion in the phase near the gap¬ 
less point is half the total multiplicity, i.e., 

M 

2d^J2mj- (C4) 

7=1 


where M is the number of zero-energy modes Ej{p*) = 0. 
Note that, due to particle-hole symmetry, M is even. 
If Ey, ^ 0, the first derivative of the Pfaffian is non¬ 
zero {d = 1) and, as a consequence, the fermion par¬ 
ity Vy = sgnEy changes at the gapless point. In 
this case Eq. (C4) and particle-hole symmetry mandates 
that there can be only two doubly-degenerate LE lev¬ 
els E±{p*) = 0 with multiplicities m± = 1, i.e., which 
have a linear phase-dispersion E±{p) oc {p — p*). This 
mandates, as we will show, a discontinuous drop in the 
Josephson CPR. A necessary condition for this case 
{Ey, ^ 0 and d = 1) is the lifting of the spin degeneracy, 
due to the broken time-reversal symmetry. If Ey, = 0 
{d > I) instead, there is no constraint in general on the 
phase-dispersion of zero-energy modes, which may be lin¬ 
ear {rrij = 1) or non-linear (m, > 1). For example, if 
d = 2, one can satisfy Eq. ( |C4| ) with M = 4 particle- 
hole and spin degenerate modes with linear dispersion 








mi = m 2 = m 3 = 7714 = 1, or with M = 2 particle-hole 
symmetric modes with parabolic dispersion mi = m 2 = 2 
if time-reversal symmetry is broken. Therefore in the 
limit points where = 0 the CPR may or 

may not be discontinuous, respectively if the zero-energy 
modes have a linear (rrij = 1) or non-linear (rrij > 1) 
phase-dispersion. 


Broken time-reversal symmetry 


In the case of broken time-reversal symmetry and if 


^ 0, one has that d = 1 in Eq. (C4) and, as a conse¬ 


quence of particle-hole degeneracy, there can be only two 
degenerate zero-energy modes E±{lp*) = 0 with multi¬ 
plicity m± = 1. Therefore one can factorize out the 


doubly-degenerate LE Andreev levels E±{^p) in Eq. (Cl) 


as 


det[HM] = |E±(^)|2 H (C5) 


where Ei{ip) are the Andreev levels with higher energy. 
Again for the properties of the determinant and of the 
Pfafhan one obtains 

E^{^) = ±^, (C6) 

Xv 

where Xv = !!£; >£;+ ^ii^) > 0 is the product of pos¬ 
itive Andreev levels with energy greater than the LE 
level. Since d^\E^\ = and sgn is 

the fermion parity, it follows that 

d^E±[ip) = ±rj^- , (C7) 

\ Xif / 


where = ^^pX^p and The contribution 

to the current of the LE level can be obtained directly 


from Eq. (C7), since Iiei<p) = {e/h)d^E-{ip). This im¬ 


plies that a change of the fermion parity at the gapless 
point mandates a sign-change of the LE level contribu¬ 
tion to the current. Moreover, at the gapless point one 


has E±[tp*) = = 0, the second term in Eq. (C7) 

vanishes and one obtains 


\d^E±{^*)\ = 




(C8) 


where Xtp* = Y\e > 0 is the product of pos¬ 

itive energy levels at tp*. The quantity Xv* 
related to the pseudo-determinant, i.e., the product of 
non-zero eigenvalues of the Hamiltonian pdet \H{ip*)] = 

considered {d = 1) one 
iias where the mi¬ 

nus sign comes from the fact that there are 2N — 1 
negative energy levels Ei{ip*) < 0. Therefore one ob¬ 
tains Xv* — \/| pdet \H{ip*)\\. In this case (d = 1, only 


two zero-energy levels at p*) the pseudo-determinant 
can be calculated using the identity p det \ 'H{ip*)\ = 
limE_,.o det + el]- Equation ( |C8| leads to 

Eq. @ by noting that the level crossing mandates a 
discontinuous drop of the Josephson current given by 
AI{ip*) = —2{e/H)\dipE±{ip*)\. Moreover, in this case 
the LE Andreev level in the neighborhood of the gapless 
point (fi* is linear in the phase tp for ^ 0, and can be 
expanded at the first order as 

E±{p)^±vJ-^{p-p*), (C9) 

Xv 

where the fermion parity Vcp changes its sign at the gap¬ 
less point p*. Therefore, discontinuities in the CPR cor¬ 
respond to fermion parity transitions in systems with bro¬ 
ken time-reversal symmetry. For the same reasons, dis¬ 
continuities of the current with respect to any parameter 
v of the Hamiltonian, e.g., magnetic field, chemical po¬ 
tential, or SOC, correspond to fermion parity transitions 
in these systems. 


Unbroken time-reversal symmetry 


In the case that time-reversal symmetry is unbroken, 
any zero-energy mode is at least four-fold degenerate due 
to particle-hole and spin degeneracy. Therefore, if E^* = 
0, there exist at least four energy levels Ej{p*) = 0 with 


multiplicity mj > 1 and, as a consequence of Eq. (C41, 


the total multiplicity is 2d = = 4n > 4. Hence, 

the lowest order d of the Pfafhan expansion in the phase 
near the gapless point is even and therefore E^, = 0. 
In this case, the CPR may still be discontinuous, but 
these discontinuities do not correspond to any topologi¬ 
cal phase transition, since the Pfafhan E^p oc {p — p*Y 
do not change its sign at the gapless point (since d is 
even). This is in agreement with the fact that an s-wave 
superconductor is topologically trivial if the time-reversal 
symmetry is unbroken. 

In the simplest case where no additional degeneracy 
is present and the multiplicity of each of the four-fold 
degenerate levels is mj = 1, one has 2 d = A and conse¬ 


quently Eq. (C5) becomes 


det [n{p)] = \E±{p)f W_ E^pf 

Ei >£/_(_ 


(CIO) 


Moreover, using again the fact that det \H{p)\ = E^ one 
obtains 


E±{p) = ±4 


\IA 

Xv 


(CIl) 


while the phase-derivative of the LE level is given by 
I ^ 


dp,E±{p) = ± 


2\/Xvl-^<, 




V E' —\E 
' '''v\ 


(CI2) 
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where = 0 at the gapless point (d = 2). Since 

the Josephson current cannot diverge, the derivative of 


the LE energy in Eq. (C12| is finite at ip*. In the limit 


of small gap A ^ t, one has = C [cos {(p — 9) + A] 
[cf. Eq. (j^J, and therefore since F^* = F'^, = 0 at the 
gapless point, one must have A = ±1. Moreover for a 
vanishing magnetic field (unbroken time-reversal symme¬ 
try) one has that 9 = mr, and therefore in this case the 
gap closes at ip* = mr for A = ±1. Finally, in the limit 
(f ^ (p*, Eq. (C12) yields 





(C13) 


Hence, in this case the CPR may still exhibit a single dis¬ 
continuity at the gapless point (p* = mr due to the pres¬ 
ence of the four-fold degenerate zero-energy mode, which 
is given by AI{(p*) = —Ae/H\d^E±{(p*)\, and which do 
not correspond to any topological phase transition. 


Appendix D: CPR at finite temperature 


At finite temperature, the Josephson current in Eq. Q 
can be rewritten using the particle-hole symmetry and 
the identity f{x) — f{—x) = — tanh(a:/ 2 ), which yield 


Hp) 


e 

h 


tanh 

Ei>0 


E^{p) 

2kBT 


d^Ei{(p). 


(DI) 


where we neglect the contribution of higher energy An¬ 
dreev levels. It follows that the phase-derivative of the 
Josephson current at low temperatures is d^pl{p*) « 
— {h/e)AI{Lp*Y/{SksT) at the gapless point, and di¬ 
verges for T —0. For this reason the scaling of the 
phase-derivative in the limit T —>■ 0 can be used to di¬ 
rectly measure the discontinuous drop AI{(p). 

We notice that the fundamental premise for the ob¬ 
servability of the effects discussed in this work is the de¬ 
tectability of the discrete nature of the subgap Andreev 
levels of the Josephson junction. This means that, in the 
ballistic limit, the average energy spacing between An¬ 
dreev levels in the low energy (subgap) spectrum should 
be larger than the other relevant energy scales of the 
system. In particular, the temperature should be suffi¬ 
ciently small, i.e., T ^ Td = Sd/ks, and disorder suffi¬ 
ciently weak, i.e., at least H/t Sd with r the scattering 
time in the nanostructure. Numerical calculations have 
confirmed that the presence of weak disorder does not 
affect the relation between the topological phase transi¬ 
tions and the discontinuities in the CPR, as long as the 
time-reversal symmetry is broken. 


Appendix E: Majorana polarization 

The total Majorana polarization[^ at any lattice site 
Xi and at energy E is defined as the complex number 


In the case of broken time-reversal symmetry, any dis¬ 
continuity in the CPR which corresponds to a topological 
phase transition is given by the contribution of a doubly- 
degenerate {d = I) Andreev level (cf. Appendix [C| . The 
contribution of the LE level E±{(p) can be considered 
separately, and one can evaluate the CPR as the sum 
of the current contributions Iie(p) and Ihe(p), respec¬ 
tively of the LE level E+{ip) > 0 and of higher energy 
levels Ei{(p) > Sd, where Sd is the energy gap between 
the lowest and the higher Andreev levels. At low tem¬ 
peratures T <^Td = Sd/kB, for any level Ei(ip) > 6d one 
has i&Tih.\Ei{ip)/{2kBT)] « 1 in Eq. (Dl), and therefore 


the contribution to the Josephson current of higher en¬ 
ergy levels does not depend on temperature, and can be 
approximated as a linear function of the phase p, i.e., 
Ihe{p) « Ihe{p*) + Ihe{p*){p “ P*) near the closing gap 
p*. On the other hand, the contribution of the LE An¬ 


dreev level can be obtained from Eq. (C9) which gives 


E+{p) « —h/{2e)AI{p*){p — p*), and therefore 


he{p) 


AI{p*) 


tanh 


h Al{p*){p-p*) 
e 4:kBT 


(D2) 


which leads to Eq. ([^. Furthermore, the first order 
phase-derivative of the current near the closing of the 
particle-hole gap yields 

H AI{p*){p - p*)' 


O T/ \ hAI{p*f 2 




(D3) 


A) = ^ 2 5{E - Ej), (El) 

j 

where the sum is over all eigenstates ipj with energy 
Ej, 5{E) the Dirac delta function, and where Uia and 
Vis are the particle-hole representation of the eigenstate 
i’j = Hi -I- + v^c\^. Note that the 

definition in Eq. ( |E1[ ) is consistent with our choice of the 
Nambu basis '5'i = (c^, C 4 , c|^^). For a pure par¬ 

ticle state {vis = 0) or pure hole (uis = 0) one has 
2{uiiv*^ + = 0 while for Majorana edge states 

localized at the lattice site in the form 

iTs = ^ , (E2) 

one has 2{uii,v*^ + Ui^v*^) = where the sign cor¬ 

responds to the two different flavors of Majorana states. 
Any linear superposition of particle and hole states with 
Uis 0 and Vis ^ 0 can be written as a superposition 
of two orthogonal Majorana states a+ 7 )^ -I- a- 7 ~. Such 
particle-hole modes have in general a finite Majorana po¬ 
larization 2{ui^v*^ + U 4 U 4 ) ^ 0 . 

At the gapless points, where the fermion parity changes 
its sign, the presence of zero-energy modes can be re¬ 
vealed by the Majorana polarization at zero energy 
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FIG. El. Majorana polarization |p(ri,0)|^ at zero energy for 
|A| < 1 in a wire as a function of the gauge-invariant phase ip 
and lattice site Vi. 


shown in Fig. El as a function of the gauge-invariant 
phase p and lattice site r^. The Majorana polarization 
is finite at gapless points inside the nanostructure. 
This entails that such zero-energy modes are not purely 
hole nor particle-like, but a superposition of orthogonal 
Majorana states. Note that these zero-energy modes are 
not localized at the edges, but are delocalized along the 
whole nanostructure. 
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